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Abstract. Solar-like oscillating giants observed by the space-borne satellites CoRoT 
and Kepler can be used as key tracers of stellar populations in the Milky Way. When 
combined with additional photometric/spectroscopic constraints, the pulsation spectra of 
solar-like oscillating giant stars not only reveal their radii, and hence distances, but also 
provide well-constrained estimates of their masses, which can be used as proxies for the 
ages of these evolved stars. In this contribution we provide supplementary material to the 
comparison we presented in Miglio et al. (2013) between populations of giants observed 
by CoRoT in the fields designated LRcOl and LRaOl. 



This contribution is organised as follows. In ^T] we briefly recall how global stellar properties 
(mass, radius, distance) are determined using average seismic constraints, while we refer to [ 1 ] for a 
general motivation behind this study, the methods adopted, and main results obtained. In ^2]we present 
the data used in this study, along with independent verifications of the asteroseismic radius and mass 



determination available in the literature {{ 2. 1 1. Comparisons between the properties of the population 
observed in LRaOl and LRcOl are presented in ^3] Finally, a comparison of observed and synthetic 
populations computed using parameterised models of the Milky Way is discussed in ^4] 



1 Determination of global stellar properties 



Radii and masses of solar-like oscillating stars can be estimated from the average seismic parameters 
that characterise their oscillation spectra: the so-called average large frequency separation {Av), and 
the frequency corresponding to the maximum observed oscillation power (v max ). We refer to [ 1 ] for a 
discussion on the method. 



a e-mail: a.miglio@bham.ac.uk 



EPJ Web of Conferences 



We determine stellar radii and masses by combining the available seismic parameters v max and Av 
with effective temperatures . The latter are determined using 2MASS [2 1 J and K s photometry and 
the colour-Teff calibrations by 0, which depend only weakly on metallicity. 2MASS colours were 
transformed into the CIT photometric system used by using the relations available in J4) and H. 

We then compute luminosities L using the Stefan-Boltzmann law: L = 4nR 2 crT^ s , where cr is the 
Stefan-Boltzmann constant. The distance modulus of each star is determined as K' M - K s q, where K s q 
is the de-reddened apparent 2MASS K s magnitude, and K' si) the absolute K s magnitude. The latter is 
obtained combining L and the bolometric corrections from [6|, which are based on Castelli & Kurucz 
2005 ATLAS9 model atmospheres, for the range of r e fj and log g into consideration. 

We take into account the effect of interstellar extinction on the magnitude and colour of each star 
using the 3D model of Galactic extinction in the V band (Ay) by Q. The extinction in the J and K s 
bandpasses were determined following |8| assuming the spectral energy distribution of a Kl giant. 
Since the extinction is distance dependent, we iterate the procedure until the derived distance does 
not vary by more than 1 %. We choose to consider magnitudes in the near-IR to reduce the effect of 
interstellar reddening in both the determination of T e g and apparent de-reddened magnitudes. 

2 Data available and uncertainties on stellar properties 

The CoRoT photometric time series used in this work were obtained in the so-called exofield during 
the first long CoRoT runs in the direction of the Galactic centre (LRcOl) and in the opposite direction 
(LRaOl). These long runs lasted approximately 140 days, providing us with a frequency resolution of 
about 0.08 //Hz. A first analysis of these data was made by J9), ifTOl and ifTTI . Solar-like oscillations 
were detected in 435 and 1626 giants belonging to LRaOl and LRcOl, respectively. To derive M,R, 
and distance we use the values (and uncertainties) of v max and Av as determined by lfT2l . The typical 
uncertainty on v max and Av for the 150-d long CoRoT observations is of the order of 2.4 % and 0.6 %, 
respectively. 

The uncertainties on the stellar properties are estimated using Monte Carlo simulations adopting 
the following constraints: 

- Apparent magnitudes: uncertainties on J and K s are taken from the 2MASS photometry available 
in the EXODAT [ 13 1 catalogue (the median uncertainty is 0.02 mag). 

- Extinction/reddening: a random error in Ay of 0.3 mag was considered both in determining r e ff 
photometrically and in de-reddening K s apparent magnitudes. 

- T e ff: for each star we considered two sources of uncertainty. The first (100 K) due to uncertainties 
on the colour- T e ff calibration itself [3|. The second one due to uncertainties on reddening. This 
results in a median combined uncertainty on T e ff of ~ 190 K (calibration+reddening). 

We notice that, when determining the distance, the effect of reddening affects not only the de- 
reddened apparent bolometric magnitude, but also the determination of L (mostly via r e ff). A higher 
A v increases the estimated T e g, hence L, but it also increases the apparent de-reddened luminosity /. 
Since d oc (L/Z) 1 ^ 2 , the overall effect of reddening on the distance itself is partly reduced (see FigfTb. 



By performing 500 realisations of the data assuming the uncertainties described above to be Gaus- 
sian, we find a median intrinsic uncertainty of 5 % in distance, 3.6 % in radius, and 10 % in mass. 

To explore the effect of possible systematic offsets in the T^g scale, we increased/decreased T e ff by 
100 K, and found that the distance estimate is affected by ~ 2.5 % (the radius by 1 % and the mass 
by 3 %). These uncertainties may be further reduced when spectroscopic constraints will be available, 
however, since the estimates of R and M are based upon seismic scaling relations, care has to be taken 
on possible systematic uncertainties which are not taken into account in the estimates presented so far. 
Scaling relations for Av and v max are based on simplifying assumptions that need to be independently 
verified, and supported with empirical tests. 

2.1 Independent verification of the asteroseismic radius/distance and mass determination 

In addition to preliminary investigations with stellar models (see e.g. lfl4l[T?j |. the Supporting Online 
Material of ifTBI . Q7]), empirical tests of the asteroseismic radius determination are now available for: 
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Fig. 1. Example of the effect of reddening on the distance determination in a LRcOl star. Left panel: change in 
r c ff due to a change in the extinction Ay. Middle panel: Absolute magnitude K' s0 as a function of Ay (blue line 
with open circles), and de-reddened apparent magnitude (red line with filled circles). The distance modulus 
determined for the central value of A v = 0.17 was subtracted from A" j0 . Right panel: Distance as a function of the 
assumed Ay. 



- 27 nearby dwarf and giant stars with precise Hipparcos parallaxes and asteroseismic constraints 
obtained with CoRoT data or ground-based spectroscopic campaigns. As shown in 1171 for this set 
of targets the mean difference between seismic radii and those estimated by classical constraints is 
-1.5%, with a standard deviation of 6 %. 

- 22 dwarfs with available Kepler asteroseismic data and precise Hipparcos parallaxes (see ifTHl ). 
The comparison shows that there is an overall agreement within + 5% between the distance ob- 
tained using the asteroseismic radii and the Hipparcos one. Additional stringent tests were carried 
out by 1 19 1, including interferometric constraints on angular radii, showing that estimates of aster- 
oseismic radii for main-sequence stars are accurate to better than 5 %. 

- Giants belonging to the old-open cluster NGC6791. There is a remarkable agreement between the 
distance determination based on eclipsing binaries [20 1 and the one based on the radius deter- 
mination of RGB stars (see I17II2U ). While the seismic and "classical' Q radii agree within 5 %, 
evidence for a small (~ 5 %) but significant discrepancy between the two determinations of radii 
was found for red-clump stars. 

- Detailed modelling of individual frequencies available so far for two RGB giants [22, 23 1 leads to 
radii and masses which agree within Icr (i.e. ~ 2% and ~ 5%, respectively) with those derived 
combining v max , Av, and T e ff. 

- The radius distribution of stars in the observed CoRoT populations is dominated by an isolated 
peak at ~ 1 1 R which is in agreement with the theoretically predicted radius of red-clump stars. 
Moreover, as reported l24l . the peak in the distribution of the absolute K s magnitudes of pulsating 
giants in the CoRoT fields is in agreement with the mean absolute K s magnitude of Hipparcos red- 
clump giants [25 26 1. This agreement represents a further argument in support of the seismically 
derived distances, since in old composite populations the absolute K s magnitude of clump stars is 
expected to depend little on the details of the population itself [see [27). 

Based on the tests reported above, we consider a 5 % systematic uncertainty on the radius determina- 
tion. This leads to an overall uncertainty on the distance of 5 % (random) + 5 % (systematic due to 
radius determination via scaling relations) + 2.5 % (systematic due to the T e ff scale). 

Independent tests of the mass determined using seismic scaling relations are less numerous and 
stringent: precise estimates of stellar mass are typically available only for visual/eclipsing binaries or 
via model-dependent estimates (e.g. isochrone fitting in clusters): 

- Direct mass measurements are available only for a handful of dwarfs with detected solar-like os- 
cillations [see 28, 29 1. Self-consistency checks of the mass determination were performed in |29| 
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combining v max or Av alone with independent radii determinations, but given the large uncertain- 
ties in the data used, these comparisons could not test the consistency of the mass determination 
(using Av or v max ) to better than 15-20 %. 

- As recalled above, [22] and [23 ] showed that, in the two cases analysed, detailed modelling of 
individual frequencies in (low-luminosity) giants leads to masses which are within 5 % those de- 
termined using scaling relations. These tests should be performed for a much wider set of targets, 
encompassing stars of different mass, metallicity, and evolutionary state. 

- The strongest constraints are provided by giants in clusters. Independent radius estimates of stars 
in NGC6791 were used in ifTTl to check consistency of the mass determination using the scaling 
relation for v max or Av alone. While no significant systematic effect was found on the RGB, IfTTll 
suggested that a relative correction to the Av scaling relation should be considered between RC 
and RGB stars, affecting the mass determination of clump stars by ~ 10 %. 

By combining constraints from near-turnoff eclipsing binaries and stellar models, [30] determined 
precisely the mass of stars on the lower RGB of NGC6719. Their value agrees within ~ 7 % with 
the average mass of RGB stars determined by GTl and IfTTll . The difference is however signif- 
icant given the quoted uncertainties, and further work is needed to understand the origin of the 
discrepancy. 

- Finally, the mass of giants in NGC68 19 derived from seismology lfT7ll2"Tl is also in good agreement 
(~ 10 %) with the one estimated from isochrone fitting OTl[32l . but the uncertainties and model- 
dependence of the latter method do not provide a stringent test to the seismically determined mass. 

It is clear that additional tests against independent estimates of mass and radius are needed, and in 
particular for metal-poor stars, which are under-represented in the above-mentioned tests. In addition 
to these empirical tests, Av computed from model frequencies may be used to correct, or at least inves- 
tigate, possible shortcomings in the Av scaling relation. Work in this direction was already presented 
in ITT31 and extended, although for a specific case only, to He-burning stars in [ 17]. This approach cap- 
tures both the deviation from the asymptotic expression due to the varying radial order of the observed 
modes, as well as to the different sound-speed distribution, hence acoustic radius, of stars with same 
mean density. In Fig.|2]we present an example of such model-predicted correction, which extend tests 
presented in iPTBI to stars in the the core-He burning phase. More detailed results will be presented in 
a forthcoming work. 

Because we have much less confidence in theoretical computations of v max - which rely on the 
complicated excitation and damping processes - than we do in theoretical predictions of the oscillation 
frequencies, model-computed v max have so far not been used in stellar properties estimation. More 
work is needed to understand the observed v max scaling, and theoretical studies have made progress on 
the problem (e.g., see 11331 ). 



3 Comparison between populations observed in LRa01 and in LRc01 

Before comparing populations of giants observed in the two fields, we check whether different biases 
were introduced in the target selection in LRaOl and LRcOl. 



3.1 CoRoT target selection 

We retrieved from EXODAT [ 1 3 1 photometric information on all the stars in the field, as well as the 
targets observed. As shown in Fig. [3] targets within each field of view were selected largely on the 
basis of colour-magnitude criteria. Within the observed targets, solar-like oscillations were searched 
for in stars belonging to a limited colour-magnitude domain: 0.6 < J - K s < 1 and K s < 12 [see 
[12). Restricting to this domain, we find no significant difference in the target selection bias applied to 
LRcOl compared to LRaOl. 
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Fig. 2. Upper panel: Ratio Av/Av B to (p/p ) 1/2 as a function of T sff for models computed using the ATON code 
[34 1, with initial heavy-elements mass fraction Z = 0.02, initial helium mass fraction Y = 0.278, and mass 
between 0.7 M Q and 3.5 M Q . Adiabatic oscillation frequencies were computed using LOSC 1351 (see 1361 1371 
for more details about the models and the computation of oscillation frequencies). (Av) is a gaussian-weighted 
average of the large frequency separation of radial modes (using a Gaussian of width v max /2 centred in v max ). 
Lower panel: HR diagram showing evolutionary tracks corresponding to the models shown in the upper panel. 



3.2 Distributions of radius and mass 

Among the 435 (1626) giants with detected solar-like oscillations in LRaOl (LRcOl), 198 (678) have 
SNR high enough not only to detect v max and Av but also to have a detailed characterisation of the 
envelope corresponding to the oscillation power excess (see Sec. 2.2 in |fj2l ). In what follows we shall 
refer to these two samples as N2 and N3, following the notation introduced in lfl2l . 

We now compare the distributions of radius and mass of stars in LRaOl and LRcOl. Since scaling 
relations were not tested at high luminosities we excluded from the sample the few stars with L > 200 
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Fig. 3. Upper-left panel: Colour-magnitude diagram of stars with R < 16 in the LRaOl field of view (blue 
circles). Stars observed by CoRoT (orange dots) in the magnitude range 0.6 < J - K s < 1 and K s < 12 (green 
dots). Upper-right panel: histogram showing the number of targets in the field (blue) and those observed (red) 
in the colour-magnitude domain selected to look for solar-like oscillating giants. The red part of the diagram 
{J — K s > 0.6) is populated mostly by giants. Lower panels: Same as in the upper panels, but for LRcOl. 



L Q , and stars with no 2MASS photometry available. The fraction of targets excluded represents < 3 % 
of the whole sample of N3 data. 

As described in |1|, we applied the Kolmogorov-Smirnov (K-S) test to the distributions of mass 
and radius to quantify differences between the two populations. We considered both N2 and N3 sets 
of data, including stars of all radii, or restricting to giants with R < 20 R . When comparing radii, the 
null hypothesis (LRaOl and LRcOl samples are drawn from the same parent distribution) cannot be 
rejected, or the difference between the two populations is marginally significant at best. On the other 
hand, we find the difference between the mass distribution of the two populations highly significant 
(K-S probability higher than 99.9 %), with stars in LRcOl having a lower average mass than those in 
LRaOl. 

To test the impact of statistical uncertainties on mass and radius on the K-S test, we perturbed the 
observed radii and masses adding a random offset drawn from a Gaussian distribution having a stan- 
dard deviation equal to the estimated uncertainty on the mass / radius. We generated 1000 realisations 
of such distributions and performed a K-S test on each realisation. The distribution of results from 
the 1000 K-S tests confirms that while comparing radius the null hypothesis cannot be rejected, the 
difference in the mass distribution is significant (in 95 % of the realisations the K-S probabilities are 
higher than 99 %) . 
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Fig. 4. TRILEGAL simulations. Distributions of radius 
(upper-left panel), mass (upper-middle panel), age (upper- 
right panel), and metallicity (lower-left panel) of giants in 
synthetic populations computed with TRILEGAL for stars in 
the field LRaOl (blue) and LRcOl (red). Stars in the syn- 
thetic population were selected in magnitude (K s < 12), colour 
(0.6 < J - K s < 1), and v max (6 /iHz < v max < 80 piiz) to ac- 
count for target selection effects. 
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Fig. 5. Same as Fig. [4] but for simulations with the Besancon Model of the Milky Way. The age sampling of the 
Besancon model is very coarse, resulting in an age resolution of only ~ 1 Gyr (see e.g. |38|) 



4 Comparison with synthetic stellar populations 

For a meaningful comparison between observed and synthetic populations, we apply to the synthetic 
population the same selection criteria based on colour and magnitude adopted in the target selection. 
Moreover, as discussed in detail in IPTZI . since no significant bias in the distribution of the targets is 
present in the v max range between 6 juHz and 80 fjHz, we shall only consider observed and simulated 
stars with v max in this frequency range. 

We consider synthetic populations computed with two codes: the Besangon Model of the Milky 
Way [see e.g.EIIEIl for a detailed description], and with TRILEGAL 01401. TRILEGAL builds a ge- 
ometric model for the Milky Way (MW), including main components such as the thin and thick disks, 
halo and bulge, each one containing a particular stellar population. The latter are defined by means of 
a star formation rate and age-metallicity relation. The main geometric parameters of the MW compo- 
nents are then calibrated using wide-area data for several lines-of-sight and for magnitudes at which 
the reddening is small and at magnitudes for which problems such as photometric incompleteness and 
star/galaxy separation are not an issue. Then, this simple model is applied everywhere in the MW, for 
different filters, and at deeper magnitudes, making the implicit assumption that the several MW com- 
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Fig. 6. TRILEGAL: stars in simu- 
lated fields centred at the galactic 
longitude of LRaOl (/ ~ 217°, left 
panel) and LRcOl (/ ~ 37°, right 
panel), and at decreasing galac- 
tic latitudes (from top to bottom). 
Stars are selected with the same 
colour-magnitude-v raax criteria as 
the observed targets (see Fig. |4j, 
and their age is colour coded. 

ponents (but for the interstellar dust) are smoothly distributed and quite uniform in their distributions 
of ages/metallicities. 

Simulations with both TRILEGAL and the Besancon Model (see Fig.|4]and[5]) show that, although 
a similar distribution of radius is expected in LRaOl and LRcOl, the age (hence mass) distribution of 
the two populations of giants is expected to differ. Stars in LRcOl are expected to be older -on average- 
than those in LRaOl. 

To check whether in the synthetic populations the difference between LRcOl and LRaOl is due to 
a different latitude, we simulated with TRILEGAL CoRoT fields at decreasing latitudes, and fixing 
the longitude to that of LRaOl and LRcOl (see Fig. [6]). The expected distribution of mass and age are 
presented in Figure [7] In the magnitude range observed by CoRoT, while a broad spectrum of age / 
mass is expected for giants in low-latitude fields, only old stars are expected to populate high-latitude 

fields. In these models the difference in the mass distribution of stars in LRcOl (b 7°) and LRaOl 

(b 2°) is dominated by the different latitudes (hence different average height below the plane). 

Our aim here is not to find the synthetic population that best matches the data, but to show that 
the difference we see in the observed distribution is in qualitative agreement with the simulations 
which include an increase of the disk scale height with age. To check the effect of this assumption 
on the predicted distributions of mass and age, we computed with TRILEGAL synthetic populations 
imposing an age-independent vertical scale height of the disk, and setting the latter to 300 pc. As 
shown in Fig.[7]in this case we do not expect to see a clear mass / age gradient with latitude. 

TRILEGAL simulations also suggest that the sample observed in the two fields is dominated by 
thin disk stars, with a combined fraction of thick disk and halo stars < 10 %. The bulge component, 
while obviously absent in LRaOl, contributes by < 1 % to the population in LRcOl. 



5 Grid based method: age estimates 

The age of RGB and red-clump stars is largely determined by their main-sequence lifetime and hence, 
to a first approximation, by their mass. To estimate the age of giants in the two observed populations 
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Fig. 7. Open bars: same as Fig. [^j showing the mass distribution for each field at different latitudes and at the 
longitude of LRaOl (left panel) and LRcOl (right panel). For reference, solid grey bars show the distribution of 
mass in simulations with no age-dependent thin-disk vertical scale height. The galactic latitude of LRcOl and 
LRaOl is b = -7° and b = -2°, respectively. Same as Fig.[6] showing the age distribution of synthetic populations 
in fields of decreasing latitudes (from top to bottom) and of the longitude of LRaOl (left panel) and LRcOl (right 
panel), with (open bars) and without (grey bars) an age-dependent vertical density scale height for the thin disk. 



we use PARAM, a Bayesian stellar parameter estimation method described in detail in [41 1. PARAM 
is largely inspired by the method developed in ll42l . which is designed to avoid statistical biases and to 
take error estimates of all observed quantities into consideration. For this work we adapted PARAM to 
include as additional observational constraints v max , Av, and, when available, the evolutionary status 
of the star (core- Helium or hydrogen-shell burning phase), as determined by the period spacing of 
gravity-dominated modes (see Il43l ). Having knowledge of the evolutionary status is particularly useful 
in the age determination of stars with an estimated radius typical of RC giants. Since the mass of core- 
He burning stars is likely to be affected by mass loss, the age-mass relation in red-clump stars may 
be different from that of RGB stars (see e.g. Fig. 1 in ifTTl ). This extended version of PARAM will be 
made available, via an interactive web form, at the URL http : // stev . oapd . inaf . it/param 

The main assumptions of the method are: a maximum value for the age, which we set to 10 Gyr, a 
prior that all ages are equally likely and, given that no information on the metallicity of these stars is 
available, a broad gaussian distribution is assumed as prior for the metallicity. The gaussian distribution 
is centred at [Fe/H]=-0.5, has a standard deviation cr = 1, and ranges from [Fe/H] = -2.3 to [Fe/H] 
= +0.3. Finally, an implicit hypothesis of the method is that theoretical models provide a reliable 
description of the way stars of different mass, metallicity, and evolutionary stage distribute along the 
red-giant phase. 

We apply this parameter estimation method to the sample of giants observed in LRcOl and LRaOl. 
The uncertainty on the age as estimated using PARAM is of the order of 30-40 %. The results obtained 
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with a model-based approach support the interpretation that the differences between the observed 
population in LRaOl and LRcOl are due to a different mass, hence age, distribution. By performing 
K-S tests on 1000 realisations of the populations taking into account Gaussian uncertainties in the 
age determination, we find the difference between the age distributions of the two populations highly 
significant. 

The precision of the age determination will significantly improve when constraints on the photo- 
spheric chemical composition will be available. In this context, and perhaps even more relevant, are 
tests to be performed on the accuracy of the age (and mass) estimates obtained assuming scaling re- 



lations for v max and Av (as discussed at the end of Sec. 2.1 1. Making use of the full set of seismic 
constraints (individual frequencies of both acoustic and gravity-dominated modes) will possibly pro- 
vide more stringent estimates of systematic uncertainties. We consider this our next, high-priority, step 
in this line of work. 
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